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Measurements of the total cross section of e^e~ —+ hadrons are presented in two different ranges of centre-of- 
mass energy. The measurements are made using the CLEO III and CLEO-c detectors at the Cornell Electron 
Storage Ring. The absolute cross sections and the values of R, the ratio of hadronic to muon pair production 
cross sections, are determined at seven centre-of-mass energies between 6.964 and 10.538 GeV. The total cross 
sections and values of R are also determined at thirteen centre-of-mass energies between 3.97 and 4.26 GeV; in 
addition, the inclusive and exclusive cross sections for D"*", D" and D'^ production are presented. Furthermore, 
for the lower centre-of-mass energy range, exclusive cross-sections are presented for final states consisting of two 
charm mesons: DD, D*D, DD* , D*W, Df D7 , 0*+ D~ , Dt D*' , DD* n and D^D'tt. 
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1. Introduction 

The determination of R, the ratio of the 
radiation-corrected hadronic cross section to the 
calculated lowest-order cross section for muon 
pair production is presented in two ranges of 
centre-of-mass energy, ^/s. The motivation for 
these studies is different in the two y/s ranges. 

The measurements of R at seven values of ^/s 
between 6.964 and 10.538 GeV test the asymp- 
totic freedom of the QCD coupling as in the range 
of y/s where u, d, s and c, quarks are produced. 
The values of a^ are extracted by comparing the 
measurements of _R to a perturbative QCD calcu- 
lation at the four- loop level [1]. 

Measurements of R in the region just above the 
cc threshold (^s = 3.970 to 4.260 GeV) exhibit a 
rich structure (for example see Ref. [5]); this re- 
flects the production of cc resonances. Interesting 
features include an enhancement at D*D* thresh- 
old (y^ — 4.02 GeV) and a broad plateau begin- 
ning at the Dl+D- threshold (V* = 4.08 GeV). 
There is considerable theoretical interest in the 
composition of these enhancements ^. However, 
there is limited experimental information avail- 
able which motivates the exclusive cross sections 
presented in this paper. The following final states 
are considered: DD, D*D, DW, D*1F, D+Dj, 
D*+Dj and D+D*', DD*tt and D*d\. The 
total cross sections and values of R are also de- 



termined at the thirteen points studies. 

The results given in this paper are presented in 
greater detail elsewhere |4l5j . 



2. CLEO III and CLEO-c 

All measurements presented are made with 
data collected at the Cornell Electron Storage 
Ring (CESR). The CLEO III detector [6] is used 
to study the data with y/s = 6.964-10.538 GeV. 
The CLEO-c detector [7] is used to study the 
data around cc threshold. The principal differ- 
ences between the CLEO III and CLEO-c detec- 
tors are the reduction in solenoid field from 1.5 
to 1.0 T and the replacement of CLEO Ill's four- 
layer silicon-strip detector by a six-layer all-stereo 
inner drift chamber. These modifications improve 
the reconstruction of low momentum charged par- 
ticles. 



Measurements 

10.538 GeV 



of i? at V^ = 6.964 



The values of \/s and integrated luminosities, 
/ Cdt, for each data point are given in Table [TJ 
Hadronic events are selected from these data by 
placing criteria on individual tracks and showers, 
as well as the whole event. 

Tracks used in the analysis are required to be 
of good quality, originate from close to the in- 
teraction point and have a momentum between 1 
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Figure 1. The upper figure compares the mea- 
surements of R presented in this paper with pre- 
vious measurements. The lower plot compares 
the results with predicted values of R; the width 
of band is given by the uncorrelated uncertainties 
on A. 



and 150% of the beam momentum. Showers must 
have at least 1% of the beam momentum and not 
be associated with a track. 

Several event variables are used to discrimi- 
nate against background. The average point of 
origin of all tracks along the beam direction is 
used to reject beam-gas, beam-wall and cosmic 
ray events. The total visible energy normalised 
to twice the beam energy and the missing mo- 
mentum in the beam direction normalised to the 
visible energy are used to remove two-photon and 
beam-gas events. Each event must contain at 
least four charged tracks and pass an event shape 
criteria to reject e+e" -^ l~^l~ {I = e,/x or r) 
events. The ratio of total calorimeter energy asso- 
ciated to tracks and isolated showers normalised 
to twice the beam energy is used to reject Bhabha 
and tau-pair events. Initial state radiation events 



Table 1 

The centre-of-mass energy, luminosity and the 
value of R for the seven measurement points with 
^/s = 6.964 - 10.538 GeV. The combined statis- 
tical and systematic uncertainty is given on the 
luminosity. The first, second and third uncertain- 
ties on R are the statistical, common systematic 
and uncorrelated systematic, respectively. 



yi (GeV) /' Cdt (pb-~ 



R 



10.538 


904.50 ± 9.00 


10.330 


149.80 ±1.60 


9.996 


432.60 ± 4.80 


9.432 


183.00 ±2.00 


8.380 


6.78 ± 0.06 


7.380 


8.48 ± 0.07 


6.964 


2.52 ± 0.02 



3.591 ± 0.003 ± 0.067 ± 0.049 
3.491 ± 0.006 ± 0.058 ± 0.055 
3.497 ± 0.004 ± 0.064 ± 0.043 
3.510 ± 0.005 ± 0.066 ± 0.037 
3.576 ± 0.024 ± 0.058 ± 0.025 
3.550 ± 0.019 ± 0.058 ± 0.020 
3.597 ± 0.033 ± 0.057 ± 0.020 



are removed by a criterion on the ratio of the 
maximum isolated shower energy normalised to 
the beam energy. 

The efficiency of the selection increases from 
82.1% at Vi = 6.964 GeV to 87.4% at ^i = 
10.538. The dominant background remain- 
ing after the event selection is from e^e~ —> 
T'^T~ events, which is estimated from simula- 
tion. Other remaining backgrounds are estimated 
to give contributions smaller than the systematic 
uncertainty assigned to the event selection. 

Radiative corrections must be applied to the 
measured total hadronic cross section to deter- 
mine R. Corrections are applied for soft photon 
radiation and vacuum polarisation. In addition, 
corrections are applied for hard initial state ra- 
diation to the continuum and lower mass qq res- 
onances. Furthermore, at three energy points, a 
correction is applied for the interference between 
a nearby T resonance and the continuum. 

Several sources of significant systematic uncer- 
tainty are considered: luminosity, radiative cor- 
rections, multiplicity corrections and event selec- 
tion criteria. The relative uncertainty on the lu- 
minosity is between 0.9% and 1.1% depending 
on ^/s. The uncertainty from the radiative cor- 
rections is dominated by those for the hadronic 
vacuum polarisation, which has an uncertainty 
of 1.0%. Some disagreement is found compar- 
ing the charged track multiplicity distribution in 



data and simulation. Therefore, the efficiency is 
determined as a function of muhiphcity before ap- 
plying to the data. This procedure is estimated to 
have an uncertainty between 0.4% and 1.4% de- 
pending on ^/s. The uncertainty related to other 
selection criteria is estimated to be between 1.0 
and 1.4% depending on y^. 

The measured values of R are given in Table [T] 
along with the associated uncertainties. The to- 
tal relative systematic uncertainty at each point 
is between 1.7 and 2.3%. The measured values 
are compared to previous measurements 18^ in the 
upper plot in Figure [1] The measurements are 
in agreement with the results from Crystal Ball, 
MDl and CLEO; however, they do not agree with 
the MARK I results. 

The value of ag is determined at each point 
using a perturbative QCD calculation of R at 
the four-loop level [T]. (This calculation ig- 
nores the quark masses; a determination includ- 
ing quark mass effects is presented in Ref. [9].) 
The compatibility of these measurements with 
others of Us is evaluated by exploiting the ex- 
pected running of as [10 , which depends on the 
QCD scale A. The measured values determine 
A = O.Sltaostoi? GeV and a(M|) = 0.126 ± 
0.005^0011' where the first uncertainties are sta- 
tistical and the second systematic. These results 
agree with the world averages [11] . The lower plot 
in Figure [T] compares the measured values of R to 
those predicted by the fitted value of A. 
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Figure 2. The sideband-subtracted momentum 
spectra for _D" (upper), Z)+ (middle) and Dg 
(lower) for data collected &i yfs = 4.17 GeV. The 
fit results for the different production mechanisms 
are also shown. 



4. Studies of exclusive charm and total 
cross sections at ^/s ~ 3.97 — 4.26 GeV 

Data collected at thirteen energy points be- 
tween ^/s = 3.97 — 4.26 GeV are studied. 
Specific criteria are used to select D"*", D^ 
and £>+ candidates; these closely follow other 
CLEO-c analyses [12]. The final states with 
the largest signal-to-background ratio are con- 
sidered for Z?+ — > K^n^TT^ and Z?° -^ K^n^. 
Eight decay modes are reconstructed to select 
D+ candidates: (j){K+ K-)tt+ , K*°{K-tt+)K+ , 
'i]{'yj)TT^ , 77(77)/3+(7r+7r*^), rj' (tt^ tt^ rj {-j j))Tr~^ , 

and Kg{TT^7T^)K'^. These modes correspond to 
16% of the total D^ decay width. 



Candidates are selected if their mass lies within 
±15 MeV of the nominal I?(s) values. Back- 
ground is subtracted using yields measured in 
I?(s) mass sidebands extrapolated into the signal 
region. The momentum distribution of the D(s) 
candidates is then used to determine the produc- 
tion mechanism. The I?°, D~^ and Df -^ (/)7r+ 
momentum distribution at ^/s — 4.17 GeV is 
shown in Figure [21 The different peaks corre- 
spond to different production mechanisms. The 
shape of distributions for the individual processes 
are determined from simulation. The data are 
fit to the different components to determine the 
yields. For Z?+ and D'^ production these are 
then corrected for efficiency, branching ratios and 
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Figure 3. The production cross sections for 

D^*^D(*') (upper), oi*'^ oi*'' (upper middle) and 
multi-body (lower middle) production. The lower 
plot shows the total charm cross section deter- 
mined by the alternate methods described in the 
text. The two-body and total production cross 
sections are compared to a model [Mj. 



luminosity to give the cross sections shown in 
the two upper plots in Figure [S] Two addi- 
tional components describing multi-body produc- 
tion D*D(*^7T are required to describe the low mo- 
mentum distribution at -^/s > 4.06 GeV. This is 
the first observation of multi-body production in 
the charm threshold region. 

Given the relative simplicity of the Dg Dg 
production and the limited statistics an alter- 
native technique is used to determine the cross 
sections. The separation of the different mech- 
anisms in the beam-energy difference (AE = 
Ejj^ — Eheam) and the beam-constrained mass 

(Mbc = V^beam-l^DA^) plaUC is UScd. Thc 



background is subtracted using AE and Mf,c side- 
bands. The resulting cross sections are shown in 
the upper middle plot in Figure [31 

The dominant sources of systematic uncer- 
tainty are the selection efficiency, yield determi- 
nation and the normalisation. Details of the un- 
certainty on the selection efficiency can be found 
in Ref. [12]. The signal functions for the deter- 
mination of D'^ and 13+ production mechanisms 
depend on the modelling of initial state radiation 
and the helicity amplitudes for D*D*; variations 
of these models over a broad range of assumptions 
leads to the systematic uncertainty. Variations in 
the Mfjc and AE selection criteria are used to esti- 
mate the uncertainty related to signal extraction 
in the D^ modes. The uncertainty on the nor- 
malisation arises from that on the measured lumi- 
nosity and the branching fractions of the modes 
reconstructed [13]. The total systematic uncer- 
tainties are between 3.4% and 6.8% for two-body 
production mechanisms; the multi-body produc- 
tion mechanisms D* Dn and D* D tt have system- 
atic uncertainties of 12% and 25% uncertainties, 
respectively. 

The results are compared to an updated calcu- 
lation of Eichten et al. [2] . There is reasonable 
qualitative agreement for most two-body produc- 
tion mechanisms apart from D* D* production in 
the yi range 4.05 to 4.20 GeV. 

The results at 4.26 GeV have the potential to 
study the nature of the y(4260). Hybrid char- 
monium [15] and tetraquark [16| interpretations 
suggest enhancements of some production mech- 
anisms; no significant enhancements are observed 
disfavouring these models. 

The sum of the exclusive cross sections should 
equal the total charm cross section. This has 
been tested with measurements using two in- 
clusive techniques. The sum of inclusive cross- 
sections for single D^, 13+ and Dg production di- 
vided by two is found to be in agreement with 
the total exclusive cross sections. In addition, 
the total hadronic cross section is determined in 
a manner similar to that described in Section [3j 
The light-quark production cross section is sub- 
tracted using measurements below cc threshold 
extrapolated with 1/s dependence. The total 
charm cross section from this method is found to 



be in agreement with the other two techniques. 
The cross sections determined by each method 
are compared to each other and a model in the 
lower plot of Figure [31 

The total hadronic cross section is radiatively 
corrected [17] to obtain the measurements of R, 
which is shown in Figured) These measurements 
are more precise and in good agreement with pre- 
vious measurements |2I18) . 
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Figure 4. Measurements of R at thirteen points 
in the cc threshold region. 



5. Conclusions 

Seven measurements of R in the range y^ = 
6.964— 10.538 GeV are used to determine a value 
of as{M'^) in agreement with the world average. 
The exclusive and inclusive charm cross sections 
are measured at thirteen ^/s values near cc thresh- 
old. Multi-body production D*D'^*'>tt is observed 
for the first time. All measurements of R are 
more precise and in agreement with earlier mea- 
surements. 
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